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T3STS OF AIRFOILS TJESIGNSD" TO DELAY THE 
COMPRESSIBILITY SUHBLE 
By John Stack 



SUIffljIARY 



Piindaraental investigations of compressibility 
phenomena for airfolla have shown that serious adverse 
changes of aerodynamic characteristics occur as the ' 
local speed over the surface exceeds the local speed of 
sound. These adverse changes have bee:n delayed to higher 
free-stream speeds by development of suitable airfoil 
shapes. The method of deriving such airfoil shapes is 
described, and aerodynamic data for a wide range of Mach 
numbers obtained from, tests of these airfoils in the 
Langley Elt.-inch high-speed tiinnel are presented. These 
airfoils, designated the NACA l6-serles, have increased 
critical Mach number. The same methods by which these - 
airfoils have been developed are amjlic^tble to other 
airplane components. ^ 

. ;. — 

INTRODUCTION .... '~ . " " 



Development of airfoil sections suitable for high- 
speed applications has generally been difficult because^ 
little was known of the: flow phenomenon that occurs at 
high speeds. A definite critical speed has been found 
at, which serious detrimental flow changes occur that 
lead to serious losses in lift and large- increases in 
drag. This flow phenomenon, called, the compressibility 
burble, was originally a propeller problem -but, with the 
development of hlgh-3pe_e_d aircraft, serious consideration 
has to be given to other parts of the. airplane. It is 
important to realize, however, that the, propeller will 
continue to offer the most serious compressibility ;,. 
problems for two reasons: first, because propeller- 
section spe.'eds-. are higiher than the speed of. the aij?plane' 
and, second.^ * because =. structural requirements .lead to 
thick sections near the root. 
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Pimdamental Inveistlgations of high-rspeed elr-flow 
phenomena recently completed (references 1 to 5) have 
provided much new inf onnation. ' Prom practical considera- 
tions an important conclusion of these investigations 
has been the determination of the critical speed, that 
is, the speed at which the corapresaibillty burble occurs. 
The critical speed was shown to be the translational 
velocity at which the sum of the translational velocity 
and the maximum local induced velocity at the surface 
of the airfoil' or other body equals the local speed of 
sound. Obviously, then, higher critical speeds can be 
attained through the development of airfoils that have 
minimum induced velocity for any given value .of the lift 
coefficient. 

Presumably, the highest critical speed will be 
attained by an airfoil that has uniform chordwlae distri- 
bution of induced velocity or, in other words, a flat 
pressure-distribution curve. All conventional airfoils 
tend to havfe high negative pressures and coi'resjjondingly ■ ■= 
high induced velocities near the "nose, which gradually 
taper off to the- air-stream conditions at the rear of 
the airfoil. ■ If the same lift coefficient can be obtained 
by decreasing the induced velocity nisar the nose and 
increasing the induced velocity over the rear portion of 
the airfoil, the critical speed Will be increased by an 
amount proportional to the decrease obtained in the 
maximum induced velocity. The ideal airfoil for' any 
given high-speed application is, then, that shape which 
at its operating lift coefficient has uniform chordwise 
distribution of induced velocity. Accordingly, an ana- 
lytical search for such airfoils has been conducted by 
raerabors of the staff of the Langley^ Memorial Aeronautical 
Laboi<atory and these airfoils have been investigated 
experimentally in the Langley 24-inch high-speed tunnel. ^ 

The first airfoils investigated" showed marked 
improvement over those shapes already available; not 
only was the critical speed increased but also the drag 
at low speeds viaa decreased considerably. Because of 
the ianrked. improvement achieved, it was considered 
desirable to extend' the thickness and the lift-coefficient 
ranges for which the original airfoils had been designed 
to obtain, data of immediate practical value before further 
extending the investigation of the fundamental aspects 
of the -problem. . . . - - 1 
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SYMBOLS . 
abscissa or .camber line 

ordinate of camber line • A- 

thickness, percent of chord 
airfoil chord 

defined by J = |(1 - cos 8 ) 
lift coefficient . 

drag coefficient " . . . " ~ - 

minimTJm drag coefficient ■ - - . . 

pitching-moment coefficient about quarter-chord 

point 

pressure coefficient 
Mach ntjmber ' . 

critical Mach number - - 

Reynolds number • '.. 

angle of attack, degrees • ' ? 

DEVELOPMENT OF AIRFOIL SERIES " • ~ 

The aerodynamic characteristics of any airfoil are',- 
in general, dependent upon the airfoil caraber line and 
the thickness form. Mean caraber lines were derived . 
analytically to obtain a unifo3?m chordwise distribution 
of Induced velocity or pressure for certain designated 
lift coefficients, and an analytical search for a thick- 
ness form that likewise has low and uniform chordwrse' 
induced-velocity distribution was then undertaken. 

Derivation of the camber line .- Glauert (referefiSS I4.) 
has derived expressions for the local induced velocity 
at a point on an airfoil (zero thickness assumed) in 
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terms of the circuletion around- an airfoil corresponding 

to a certain distribution of vorticity along the airfoil 
surface. By assuming the distribution of vorticity to be 
constant, a line airfoil is determined that gives unlforoi 
chordwise pressure distribution. The form of the equa-i- 
tion ao derived is 

Is. = ^ Clog -X-, + X log l-r^LJL) (1) 



where y^ is the ordinate of the mean camber line, x is 
the abscissa, and the chord is taken as unity. The 
idealized form described by this equation has disconti- 
nuities at the nose and at the tail:, .This difficulty la. 
circumvented by assuming very slight gradients in the 
chord>vise load distribution just at the nose and just at. 
the tail. This forin, derived by ysing the Fourier series 
method, is given by the equation 

^= H^(0.3835 - 0.3555 COS 29 - 0.0333 cos kd. 

- 0.0095 cos 6e - 0 ,00/10 cos 89- 

- 0.0020 cos 109 - 0.0012 ^30^126) (2) 

where — = i-(l - oos 0) and c is the airfoil chord. 

Equation (2) expresses the mean cambist line choSBn for 
airfoils of the series developed. Load or induced- 
velocity ^y^adings derived from both equations (1) and (2) 
are actually fdentical for all practical . purposes . Mean- 
camber-line ordinates are given in table I for C-^ - 1.0. 

In- or-der to . obtain tlas mean camber ,line .giving uniform 
chordwise , distrifeutlon 'of- inducted velocity for other 
values of the -lif t cpefi'icient , the values glvsn.'iri . 
table I are multiplied by the value of the -desired .lift 
coefficient. ' 

E e-rivation ,of the." thickness form .- The derivation 
of the .thickness, form ia -not §q, simple or direct.as the 
derivation of the, mean camber line. The theoretical 
pressure^-dis trlbufcioa,w.as;,computed by the methods of 
reference-. 5 ^.f or, each;-Qirj th"e~:."s-6veral thickness- forms 
investigated.- in reference .6.^ . .Some of these .forms 
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approached the desired shape but further modifications 
were investigated analytically and, finally, two shapes 
were chosen for tests. These shapes, the HACA O7-OO9 
and the NACA I6-OO9, and the thporetic?-! pressure distrl-- 
bution for each are shown In f igur.e 1. .The complete 
airfoil profile is derived by first cAlculating the mean 
camber line for the desired lift coefficient and then 
laying out the thickness ordinates given in table II from 
the camber line along perpendiculars to this line. 

Airfoil designation .- Because the ideal series of 
airfoils requires an extremely large variation of shape, 
it becomes practically impossibler to use previous nvaa- 
bering systems and, further, because this new series of 
airfoils is designed to obtain a specific pressure diagram, 
these airfoils are designated by a new series of 'numbers 
that is related to the flow and the operating character- 
istics of the airfoil. The first number is a serial 
number that describes the class of pressure distribution, 
the second nvimber gives the location of the maximum nega- 
tive pressure in percent of chord from the leading edge,, 
the first number following the dash gives Hie lift coef-' 
ficlent for which the airfoil was designed to. operate, ' ' " 
and the last two hxambers give the airfoil thickness in 
"percent of chord. Thus the NACA I6-509 airfoil has "the 
shape of the NACA I6-OO9 disposed about the uniform chord- 
wise load camber line designed for a lift coSfficient of 

0.5* -- - 

Airfoils investigated ,- As previous-ly stated, two 
basi5 airfoils were investigated. (See fig. 1.) The 
NAPA 07-009 airfoil should, theoretically, give higher 
critical speed than the ITACA I6-OO9 but an earlier inves- 
tigation (reference 3)- indicated that, for pressures 
occurring near the leading edge, the" increase in the 
pressure coefficient as a result of compressibility 
effects was greater than that for pressures occurring 
farther back on the airfoil. Consequently, it was believed 
that, at speeds as high as the critical speed, the NACA 
07-009 airfoil "might have, as a result of compressibility 
effects, a pressure peak near the leading edge." " "fhe' 
NACA 16-009 airfoil was therefore developed in an attempt, 
to achieve the- uniform chordwise load distribution at 
high speeds. Both forms were tested arid the results 
showed higher drag and lower critical speed for the NACA 
07-009 airfoil. Accordingly, the NACA I6-OO9 airfoil 
was chosen as the basic form for a aeries of airfoils 
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designed to operate, at various lift coefficients. . For 
one value of the lift coefficient the effect of thickness 
variation was also investigated. The alrf.olls tested,- 
of which profiles are shown in figiire 2," are as follows: 

NACA 16-009 NACA l6-50^ 
NACA 16-109 - NACA I6-512 ^ 

NACA 16-209 NACA I6-515 

NACA 16-509 NACA I6-52I 

NACA 16-709 NACA I6-53O 

NACA 16-1009 ■ NACA I6-IO6 

NACA 07-009 NACA O7-509 



APPARATUS AND METHOD 



The teats were conducted in the Langley 2i<.-inch . 
high-speed tunnel, in which velocities _at>proaching the 
speed of sotmd can \)e. pl?tained. . A brjef description of 
this tunnel is given in reference 3« The balance meas- 
urei3 lift, drag, and pitching qjoment £nd, except for 
improvements .that permit a more accurate determlnetlon 
of the forces, is similar in principle to the balance 
used in the Langley llr.inch high-speed tunnel. The 
methods of operation are likewise similar to those 
employed in the operation of - the Langley 11-inch high- 
speed tunnel (reference 7)' 

The models were of 5 -inch chord and JO-inch span 
and were made' of • dur&lumin. A complete description of 
the nethod of constructing the models la given in ref- 
erence 8. The model mounting is similar to that used 
in the -Langley 11-inch high-speed tunnel (reference "J). 
The ttodel extends across the tunnel and through holes, 
which are of the same shape as but slightly larger than 
the model, cut in flexible brass end plates that preserve 
the contour of the tunnel walls. The model ends are 
securisd in the balance, which extends halfway around the 
test section and is enclosed in the airtight tunnel 
chamber similar to the installation in the Langley 11-lnch 
high-speed tunnel (reference J), ... 

The speed range oyer which keasuT-ements were made 
extended, in general, from 25 percent of the speed of 
sound to values in excess of th-e critical speed. The 
corresponding Reynolds number range waa i*rom approxi- 
mately 700,000 to nearly 2,000,000. The lift-coefficient 
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range for which tests were made ex;tended from zero lift 
for each airfoil to values approaching maximum lift. 



Accidental 'errors are indicated by the scatter on 
the plots showing the measured test data (figs. 5 and 
25 to 35). These errors are, in general, rather siiisll 
and affect neither the application nor the comparison of 
the data. Tunnel effects arising from end lealqage, 
restriction, and the usual type of tunriel-wall 'effect 
are important. Exact knowledge of these various effects 
is incomplete at the present tim's. The largest 'effects 
appear to arise from air leakage through the clearance 
between the model ^nd the brass end plates in the tunnel 
wall through which the model passes. Investigations of 
the leakage effects have been made for the NACA 0012 air- 
foil with a special .type of internal gap or clearance - 
that permits wide variation of the gap. . Data "obtained' 
with various gap Settings of 0.01 inch and larger Bxtra- 
polate'd to zero gap were used to evaluate the leakage 
correction for the standard type of aoixnting.. These 
corrected data were then checked by means of- wake-survey 
drag measurements with end leakage eliminated by rubbed" 
seals. Because the balance chamber is Sirtight, the 
end-leakage condition is related to. the presBUrs distri- 
bution around the' model. It was therefore considered 
advisable to check .the method of corr'ection for' end/ 
leakage by" wake-survey tests with end leaka'gfe eliminated' ' 
by rubber seals for these new airfoils, which 'have rsrdt- 
cally different pressure distributions from the older 
airfoils such as the NACA 0012. Some of these data are 
shown in figure 5. In general," the -agreement is excel- 
lent. The data have accordingly been corrected for end- 
leakage effects. ■ . ■ 

Other tunnel effects have not been completely inves- 
tigated and the data have not been corrected for such 
effects as restriction or the' more usual type of wall 
effects As presented, the data are therefore conseirva- ' 
tive, .inasmuch as investigations made thus far' indicate 
that the coefficients are high and the critical speeds 
may be low. Strictly comparable data fojp two older air- 
foils,- the 5C8 and the NAGA 2kO^-3ki for'/two^ Mach ikimbiers 
are included so that-comparlsgns can be- made. — - 
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Aerodynamic charscteristlcs for several or the NACA 
16 -series airfolXs are given in figures i{. to I8. Exajai- 
nation of these figures indicates two important discrep- 
ancies between the theoretical design conditions and the 
data obtained from the tests: First, none of the' airfoils 
attains the design lift coefficient at the desijjn angle 
of attack (0°) and/ second, the departure increa'i^es 
m&i'i^edlj with the design lift coefficiant. The departures 
radf_ be important- if varlstio'n. from the ideal pressure 
distribution is rapid- with change in lift coefficient. 
This effect, if great, would tend to cause lower drag ; 
apd higher critical speed far a narrow region near the- 
dfe sign-, condition than are shown by these data. These ' 
departures also increase with the airfoil thickness. 

The .differences between the design- coiidltions and 
^the actual test results may be expected because of the ' 
"s'lmplif ying assumptions of the thin-airfoil theory. 
Theoretically, it -Is assumed that the induced -velocities 
are negligibly small as compared with the stream velocity. 
For thin airfoils at low lifts, this approximation is 
valid. With increases of Lift or thickness, however, the 
induced velocltias approach and sometimes .exceed the 
stre-am velocity. Study of these effects, appears to-be 
very important in order to obtain the proper airfoils" 
for high lift coefficients and large thickness ratios. 
Deviations shown by the airfoils in this series having 
high lift and high thickness ratios appear to indicate 
that, the use of "a single basic shape is. ..unv/arr anted if 
it 1,3^ 'desired to obtain optimum ;airf oils .for a wide range 
oi'' I'i'fjt; .Cipeff icient and thicknesjs ./distribution.- 

Theoretical pressur^-distrij)utic»n. diagrams for the 

thlck<?r airfoils 'Showed much greater slope, of the pres- 
sure "curve than is shov/n by the basic NACA I6-OO9 airfoil. 
Preliminary study indicatied that. Incr.easirsg the leading- 
edge radius -aAfi-.-the fullness of th'e .'airfoil between the 
leading edge .and :.the maximum ordinate .niay lead to con- 
sider abrle. iniproveirientr-T5V"SJ? the thicker airfoils hererin 
reporte'^d, - • ■ 

0 omp ai? 1 s.oh of ulff o 1 1 s ^ - Plg.ure3 I9 and 20 .illua trate 
the. differences in aerodynamic characteristics between 
older propeller-blsde sections and the NACA l6-aeries 
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airfoils. At lower ■ speeds (M - Q'k5 >■ fig- 19) the 
5C8 airfoil appears to attain a much higher fflaximuin lift 
coefficient than the new airfoils. This result is impor- 
tant in that the wider useful angle-of -attack range may 
frequently be required to prevent stalling af a propeller 
during take-off. Over the normal flight range,- however, 
and in most cases for which rational choice of section 
can be made, the lower drag of the "new sections offers 
considerable opportunity to achieve higher efficiencies. 
The low drag attained by the NACA 2ii.09-3ti. airfoil developed 
from earlier tests in the Langley 11-inch high-speed tun- 
nel may appear surprising. Actually the type of flow for 
this airfoil approaches the flow that might be expected 
for the NACA I6-309 airfoil. - • . ■ 

The low drag common to most, of the NACA l6-series 
airfoils is associated with more extensive regions of 
laminar flow in the boundary layer resulting from the 
rearward position of the point of maximrum negative pres- 
sure. Unf ortxinately, however, the Reynolds number is so 
low that effects of laminar separation may appear and 
some pressure drag might occur. The small differences 
in drag between the, envelope polar for the new airfoils 
and for the NACA 2l+69-5ij. airfoil are probably a result 
of this phenomenon. Actually tiie point of_ maximum nega- 
tive pressure for the new airfoils is considerably 
farther back than the corresponding point for the NACA 
2ii.09-5i4. airfoil but, if laminar separation occurs early, 
nearly equal drag coefficients might be expected. 

At high speeds (M = O.75, fig. 20), the region Tor 
which the NACA l6-series airfoils were" desig^ied, the 
superiority of the new airfoils is clea'r.' ~The"' "earlier 
onset of the compressibility ef fects ,f or the older air- 
foils leads to early drag increases and lowered maximum 
lift coefficients. At speeds above M = 0.75 the use of 
the 6lder sections appears unwarranted for any purpose. 

Critical speed .- The variation of the critical 
speed with lift coefficient and with thickness is given 
in figures 21 and 22, respectively. " These curves indicate 
that critical speeds exceeding the theoretical values 
were attained in the tests. In the choice of the test 
critical speeds, the values were selected on the basis 
of earlier experience that indicated some rise in drag 
before large flow disturbances occurred. If these speeds- 
were chosen as the highest values reached before any 
appreciable drag; increment occurred, the agreement with 
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the theoretical .curves would be very gao'd. For ooniparlson 
the critical speed of the 3c8 aii-roll is plotted. in £a.g-- 
ure 21". The difference between .the new and the older 
airfoils is greater than shown by the curves because the 
5C3 is 8 percent thick, or 1 percent of 'the chord thinner 
than the airfoils pf the NACA I6-609 series. 

Minimum drag . - Coefficients of minimum drag plotted 
ageilnst Reynolds number are given in figures ,25 and 2J4.. 
^he lowest drag coefficient was obtained for the NACA 
16-106 airfoil; this coefficient is approximately 0,0026 
at low speeds and increases ttr-approximately O.OO32 
immediately below the critical gpeed. Cf the 9-percent- 
thick airfoil series designed to operate at various lift 
coefficients^ i?he NACA I-6-IO9 airfoil appears to have 
the lowest drag. This result is contrary to expectation 
because the symmetrical or basic form of the NACA 16-OO9 
would normally have the lov/est minimum drag ooeff iciient. 
The "difference may be due to some irregularity of the 
airfoil surface. 

The comparison of the. minimum drag coefficients 
for the 3C8 and the NACA 2i|.09-3i^. and 16-209 airfoils is 
shown in figure 25 . The high critical speed -for the NACA 
16-209 airfoil is apparent.. The comparison as given- 
directly by figure 25 is a little misleading because of 
the smaller thickness ratio for the JCS airfoil. For 
equal thickness ratios, the differences between the 
C-series and the NACA l6-series airfoils will be greater 
than shown. 

Use of the data . - The envelope polars that may be 
drawn for the NACA I6 -series airfoils represent a "new 
and much lower drag as well as higher critical speed "■ 
attainable for the design of propeller-blade sections. ' 
Even though the angle-pf -attack range Is less than for 
the older sections, there will be numerous designs for' 
which sufficient angle-of -attack range is given by the 
new sections, ^ov high-speed, high-altitude aircraft, 
the advantages of the low drag and hig^i critical speed 
are of paramount Importance and, in these designs, rational 
choice of section is of increasing" Importance . In many 
designs the diameter is fixed by considerations other ■ 
than p3r»opeller efficiency. Thus the induced losses are 
fixed and propellers of hlgjiest efficiency can be developed 
only by operating and designing- the blade sections to 
operate on the envelope polfsr^s. Another important con- 
sideration In using new blade sections to achieve highest 
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efficiency concerns the adaptation of the sections to 
older propeller designs. Optimum efficiency cannot be 
achieved by simply substituting the new sections for the 
old on a givBn design. The use of better blade . sfeetions 
permits the use of larger diameter and necessitates some 
plan-form changes., .All these factors should be .considered 
in a design for bes't" eTf iciency with the new blade" sections. 



CONCLUSION 



By a new approach to airfoil design based upon 
findings of fundamental flow studies, a new series of 
airfoils, the NACA I6 series, have been develoged which 
have increased critical Mach number and at low speeds 
reduced drag. ... .. . - ,"- _ Jr-:~-^■^-.•-~^- 
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TABI^I I 



CAMHSR-LINE ORDINATES FOR NACA I6- AND 07-S3RIES 
AIRFOILS V/ITH = -L.p ^ 

[All values measured in percent chord 
from chord Ixn^ 



Station 



0 

1-25 

^ 5. Oft 

7.5^ 
,10.00 

.15.00 

20^00 

.25.00 

; ^Q.oo 
14.0-00 
50-00 

60 -00 
70 .00 

,90.00 
95-06 
100.60 



Ordinate 



0 

.930 
1.580 

2.120 
2. 587 

5.556 

5,516 

4,861 
5.982 

2.5G7 
1.5S0 



Slope 



^3gT. 
.2^155 

.lSfe95 

,m66 

.15601+ 
. 110^2 

.087U3 

.0671+5 
.03227 

0 

- .03227 

- .067i+5 
^ .1105 2 

- .17W 

- .25i|-52 

- .6225I4. 



NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 



NATIONAL ADVISORY COmiTTSS FOR AERONAtJTlCS"- 



TECHNICAL NOTE HQ, 976 



T3STS OF AIRFOILS TJESIGNSD" TO DELAY THE 
COMPRESSIBILITY SUHBLE 
By John Stack 



SUIffljIARY 



Piindaraental investigations of compressibility 
phenomena for airfolla have shown that serious adverse 
changes of aerodynamic characteristics occur as the ' 
local speed over the surface exceeds the local speed of 
sound. These adverse changes have bee:n delayed to higher 
free-stream speeds by development of suitable airfoil 
shapes. The method of deriving such airfoil shapes is 
described, and aerodynamic data for a wide range of Mach 
numbers obtained from, tests of these airfoils in the 
Langley Elt.-inch high-speed tiinnel are presented. These 
airfoils, designated the NACA l6-serles, have increased 
critical Mach number. The same methods by which these - 
airfoils have been developed are amjlic^tble to other 
airplane components. ^ 

. ;. — 

INTRODUCTION .... '~ . " " 



Development of airfoil sections suitable for high- 
speed applications has generally been difficult because^ 
little was known of the: flow phenomenon that occurs at 
high speeds. A definite critical speed has been found 
at, which serious detrimental flow changes occur that 
lead to serious losses in lift and large- increases in 
drag. This flow phenomenon, called, the compressibility 
burble, was originally a propeller problem -but, with the 
development of hlgh-3pe_e_d aircraft, serious consideration 
has to be given to other parts of the. airplane. It is 
important to realize, however, that the, propeller will 
continue to offer the most serious compressibility ;,. 
problems for two reasons: first, because propeller- 
section spe.'eds-. are higiher than the speed of. the aij?plane' 
and, second.^ * because =. structural requirements .lead to 
thick sections near the root. 
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Pimdamental Inveistlgations of high-rspeed elr-flow 
phenomena recently completed (references 1 to 5) have 
provided much new inf onnation. ' Prom practical considera- 
tions an important conclusion of these investigations 
has been the determination of the critical speed, that 
is, the speed at which the corapresaibillty burble occurs. 
The critical speed was shown to be the translational 
velocity at which the sum of the translational velocity 
and the maximum local induced velocity at the surface 
of the airfoil' or other body equals the local speed of 
sound. Obviously, then, higher critical speeds can be 
attained through the development of airfoils that have 
minimum induced velocity for any given value .of the lift 
coefficient. 

Presumably, the highest critical speed will be 
attained by an airfoil that has uniform chordwlae distri- 
bution of induced velocity or, in other words, a flat 
pressure-distribution curve. All conventional airfoils 
tend to havfe high negative pressures and coi'resjjondingly ■ ■= 
high induced velocities near the "nose, which gradually 
taper off to the- air-stream conditions at the rear of 
the airfoil. ■ If the same lift coefficient can be obtained 
by decreasing the induced velocity nisar the nose and 
increasing the induced velocity over the rear portion of 
the airfoil, the critical speed Will be increased by an 
amount proportional to the decrease obtained in the 
maximum induced velocity. The ideal airfoil for' any 
given high-speed application is, then, that shape which 
at its operating lift coefficient has uniform chordwise 
distribution of induced velocity. Accordingly, an ana- 
lytical search for such airfoils has been conducted by 
raerabors of the staff of the Langley^ Memorial Aeronautical 
Laboi<atory and these airfoils have been investigated 
experimentally in the Langley 24-inch high-speed tunnel. ^ 

The first airfoils investigated" showed marked 
improvement over those shapes already available; not 
only was the critical speed increased but also the drag 
at low speeds viaa decreased considerably. Because of 
the ianrked. improvement achieved, it was considered 
desirable to extend' the thickness and the lift-coefficient 
ranges for which the original airfoils had been designed 
to obtain, data of immediate practical value before further 
extending the investigation of the fundamental aspects 
of the -problem. . . . - - 1 
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SYMBOLS . 
abscissa or .camber line 

ordinate of camber line • A- 

thickness, percent of chord 
airfoil chord 

defined by J = |(1 - cos 8 ) 
lift coefficient . 

drag coefficient " . . . " ~ - 

minimTJm drag coefficient ■ - - . . 

pitching-moment coefficient about quarter-chord 

point 

pressure coefficient 
Mach ntjmber ' . 

critical Mach number - - 

Reynolds number • '.. 

angle of attack, degrees • ' ? 

DEVELOPMENT OF AIRFOIL SERIES " • ~ 

The aerodynamic characteristics of any airfoil are',- 
in general, dependent upon the airfoil caraber line and 
the thickness form. Mean caraber lines were derived . 
analytically to obtain a unifo3?m chordwise distribution 
of Induced velocity or pressure for certain designated 
lift coefficients, and an analytical search for a thick- 
ness form that likewise has low and uniform chordwrse' 
induced-velocity distribution was then undertaken. 

Derivation of the camber line .- Glauert (referefiSS I4.) 
has derived expressions for the local induced velocity 
at a point on an airfoil (zero thickness assumed) in 
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terms of the circuletion around- an airfoil corresponding 

to a certain distribution of vorticity along the airfoil 
surface. By assuming the distribution of vorticity to be 
constant, a line airfoil is determined that gives unlforoi 
chordwise pressure distribution. The form of the equa-i- 
tion ao derived is 

Is. = ^ Clog -X-, + X log l-r^LJL) (1) 



where y^ is the ordinate of the mean camber line, x is 
the abscissa, and the chord is taken as unity. The 
idealized form described by this equation has disconti- 
nuities at the nose and at the tail:, .This difficulty la. 
circumvented by assuming very slight gradients in the 
chord>vise load distribution just at the nose and just at. 
the tail. This forin, derived by ysing the Fourier series 
method, is given by the equation 

^= H^(0.3835 - 0.3555 COS 29 - 0.0333 cos kd. 

- 0.0095 cos 6e - 0 ,00/10 cos 89- 

- 0.0020 cos 109 - 0.0012 ^30^126) (2) 

where — = i-(l - oos 0) and c is the airfoil chord. 

Equation (2) expresses the mean cambist line choSBn for 
airfoils of the series developed. Load or induced- 
velocity ^y^adings derived from both equations (1) and (2) 
are actually fdentical for all practical . purposes . Mean- 
camber-line ordinates are given in table I for C-^ - 1.0. 

In- or-der to . obtain tlas mean camber ,line .giving uniform 
chordwise , distrifeutlon 'of- inducted velocity for other 
values of the -lif t cpefi'icient , the values glvsn.'iri . 
table I are multiplied by the value of the -desired .lift 
coefficient. ' 

E e-rivation ,of the." thickness form .- The derivation 
of the .thickness, form ia -not §q, simple or direct.as the 
derivation of the, mean camber line. The theoretical 
pressure^-dis trlbufcioa,w.as;,computed by the methods of 
reference-. 5 ^.f or, each;-Qirj th"e~:."s-6veral thickness- forms 
investigated.- in reference .6.^ . .Some of these .forms 
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approached the desired shape but further modifications 
were investigated analytically and, finally, two shapes 
were chosen for tests. These shapes, the HACA O7-OO9 
and the NACA I6-OO9, and the thporetic?-! pressure distrl-- 
bution for each are shown In f igur.e 1. .The complete 
airfoil profile is derived by first cAlculating the mean 
camber line for the desired lift coefficient and then 
laying out the thickness ordinates given in table II from 
the camber line along perpendiculars to this line. 

Airfoil designation .- Because the ideal series of 
airfoils requires an extremely large variation of shape, 
it becomes practically impossibler to use previous nvaa- 
bering systems and, further, because this new series of 
airfoils is designed to obtain a specific pressure diagram, 
these airfoils are designated by a new series of 'numbers 
that is related to the flow and the operating character- 
istics of the airfoil. The first number is a serial 
number that describes the class of pressure distribution, 
the second nvimber gives the location of the maximum nega- 
tive pressure in percent of chord from the leading edge,, 
the first number following the dash gives Hie lift coef-' 
ficlent for which the airfoil was designed to. operate, ' ' " 
and the last two hxambers give the airfoil thickness in 
"percent of chord. Thus the NACA I6-509 airfoil has "the 
shape of the NACA I6-OO9 disposed about the uniform chord- 
wise load camber line designed for a lift coSfficient of 

0.5* -- - 

Airfoils investigated ,- As previous-ly stated, two 
basi5 airfoils were investigated. (See fig. 1.) The 
NAPA 07-009 airfoil should, theoretically, give higher 
critical speed than the ITACA I6-OO9 but an earlier inves- 
tigation (reference 3)- indicated that, for pressures 
occurring near the leading edge, the" increase in the 
pressure coefficient as a result of compressibility 
effects was greater than that for pressures occurring 
farther back on the airfoil. Consequently, it was believed 
that, at speeds as high as the critical speed, the NACA 
07-009 airfoil "might have, as a result of compressibility 
effects, a pressure peak near the leading edge." " "fhe' 
NACA 16-009 airfoil was therefore developed in an attempt, 
to achieve the- uniform chordwise load distribution at 
high speeds. Both forms were tested arid the results 
showed higher drag and lower critical speed for the NACA 
07-009 airfoil. Accordingly, the NACA I6-OO9 airfoil 
was chosen as the basic form for a aeries of airfoils 
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designed to operate, at various lift coefficients. . For 
one value of the lift coefficient the effect of thickness 
variation was also investigated. The alrf.olls tested,- 
of which profiles are shown in figiire 2," are as follows: 

NACA 16-009 NACA l6-50^ 
NACA 16-109 - NACA I6-512 ^ 

NACA 16-209 NACA I6-515 

NACA 16-509 NACA I6-52I 

NACA 16-709 NACA I6-53O 

NACA 16-1009 ■ NACA I6-IO6 

NACA 07-009 NACA O7-509 



APPARATUS AND METHOD 



The teats were conducted in the Langley 2i<.-inch . 
high-speed tunnel, in which velocities _at>proaching the 
speed of sotmd can \)e. pl?tained. . A brjef description of 
this tunnel is given in reference 3« The balance meas- 
urei3 lift, drag, and pitching qjoment £nd, except for 
improvements .that permit a more accurate determlnetlon 
of the forces, is similar in principle to the balance 
used in the Langley llr.inch high-speed tunnel. The 
methods of operation are likewise similar to those 
employed in the operation of - the Langley 11-inch high- 
speed tunnel (reference 7)' 

The models were of 5 -inch chord and JO-inch span 
and were made' of • dur&lumin. A complete description of 
the nethod of constructing the models la given in ref- 
erence 8. The model mounting is similar to that used 
in the -Langley 11-inch high-speed tunnel (reference "J). 
The ttodel extends across the tunnel and through holes, 
which are of the same shape as but slightly larger than 
the model, cut in flexible brass end plates that preserve 
the contour of the tunnel walls. The model ends are 
securisd in the balance, which extends halfway around the 
test section and is enclosed in the airtight tunnel 
chamber similar to the installation in the Langley 11-lnch 
high-speed tunnel (reference J), ... 

The speed range oyer which keasuT-ements were made 
extended, in general, from 25 percent of the speed of 
sound to values in excess of th-e critical speed. The 
corresponding Reynolds number range waa i*rom approxi- 
mately 700,000 to nearly 2,000,000. The lift-coefficient 
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range for which tests were made ex;tended from zero lift 
for each airfoil to values approaching maximum lift. 



Accidental 'errors are indicated by the scatter on 
the plots showing the measured test data (figs. 5 and 
25 to 35). These errors are, in general, rather siiisll 
and affect neither the application nor the comparison of 
the data. Tunnel effects arising from end lealqage, 
restriction, and the usual type of tunriel-wall 'effect 
are important. Exact knowledge of these various effects 
is incomplete at the present tim's. The largest 'effects 
appear to arise from air leakage through the clearance 
between the model ^nd the brass end plates in the tunnel 
wall through which the model passes. Investigations of 
the leakage effects have been made for the NACA 0012 air- 
foil with a special .type of internal gap or clearance - 
that permits wide variation of the gap. . Data "obtained' 
with various gap Settings of 0.01 inch and larger Bxtra- 
polate'd to zero gap were used to evaluate the leakage 
correction for the standard type of aoixnting.. These 
corrected data were then checked by means of- wake-survey 
drag measurements with end leakage eliminated by rubbed" 
seals. Because the balance chamber is Sirtight, the 
end-leakage condition is related to. the presBUrs distri- 
bution around the' model. It was therefore considered 
advisable to check .the method of corr'ection for' end/ 
leakage by" wake-survey tests with end leaka'gfe eliminated' ' 
by rubber seals for these new airfoils, which 'have rsrdt- 
cally different pressure distributions from the older 
airfoils such as the NACA 0012. Some of these data are 
shown in figure 5. In general," the -agreement is excel- 
lent. The data have accordingly been corrected for end- 
leakage effects. ■ . ■ 

Other tunnel effects have not been completely inves- 
tigated and the data have not been corrected for such 
effects as restriction or the' more usual type of wall 
effects As presented, the data are therefore conseirva- ' 
tive, .inasmuch as investigations made thus far' indicate 
that the coefficients are high and the critical speeds 
may be low. Strictly comparable data fojp two older air- 
foils,- the 5C8 and the NAGA 2kO^-3ki for'/two^ Mach ikimbiers 
are included so that-comparlsgns can be- made. — - 
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Aerodynamic charscteristlcs for several or the NACA 
16 -series airfolXs are given in figures i{. to I8. Exajai- 
nation of these figures indicates two important discrep- 
ancies between the theoretical design conditions and the 
data obtained from the tests: First, none of the' airfoils 
attains the design lift coefficient at the desijjn angle 
of attack (0°) and/ second, the departure increa'i^es 
m&i'i^edlj with the design lift coefficiant. The departures 
radf_ be important- if varlstio'n. from the ideal pressure 
distribution is rapid- with change in lift coefficient. 
This effect, if great, would tend to cause lower drag ; 
apd higher critical speed far a narrow region near the- 
dfe sign-, condition than are shown by these data. These ' 
departures also increase with the airfoil thickness. 

The .differences between the design- coiidltions and 
^the actual test results may be expected because of the ' 
"s'lmplif ying assumptions of the thin-airfoil theory. 
Theoretically, it -Is assumed that the induced -velocities 
are negligibly small as compared with the stream velocity. 
For thin airfoils at low lifts, this approximation is 
valid. With increases of Lift or thickness, however, the 
induced velocltias approach and sometimes .exceed the 
stre-am velocity. Study of these effects, appears to-be 
very important in order to obtain the proper airfoils" 
for high lift coefficients and large thickness ratios. 
Deviations shown by the airfoils in this series having 
high lift and high thickness ratios appear to indicate 
that, the use of "a single basic shape is. ..unv/arr anted if 
it 1,3^ 'desired to obtain optimum ;airf oils .for a wide range 
oi'' I'i'fjt; .Cipeff icient and thicknesjs ./distribution.- 

Theoretical pressur^-distrij)utic»n. diagrams for the 

thlck<?r airfoils 'Showed much greater slope, of the pres- 
sure "curve than is shov/n by the basic NACA I6-OO9 airfoil. 
Preliminary study indicatied that. Incr.easirsg the leading- 
edge radius -aAfi-.-the fullness of th'e .'airfoil between the 
leading edge .and :.the maximum ordinate .niay lead to con- 
sider abrle. iniproveirientr-T5V"SJ? the thicker airfoils hererin 
reporte'^d, - • ■ 

0 omp ai? 1 s.oh of ulff o 1 1 s ^ - Plg.ure3 I9 and 20 .illua trate 
the. differences in aerodynamic characteristics between 
older propeller-blsde sections and the NACA l6-aeries 
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airfoils. At lower ■ speeds (M - Q'k5 >■ fig- 19) the 
5C8 airfoil appears to attain a much higher fflaximuin lift 
coefficient than the new airfoils. This result is impor- 
tant in that the wider useful angle-of -attack range may 
frequently be required to prevent stalling af a propeller 
during take-off. Over the normal flight range,- however, 
and in most cases for which rational choice of section 
can be made, the lower drag of the "new sections offers 
considerable opportunity to achieve higher efficiencies. 
The low drag attained by the NACA 2ii.09-3ti. airfoil developed 
from earlier tests in the Langley 11-inch high-speed tun- 
nel may appear surprising. Actually the type of flow for 
this airfoil approaches the flow that might be expected 
for the NACA I6-309 airfoil. - • . ■ 

The low drag common to most, of the NACA l6-series 
airfoils is associated with more extensive regions of 
laminar flow in the boundary layer resulting from the 
rearward position of the point of maximrum negative pres- 
sure. Unf ortxinately, however, the Reynolds number is so 
low that effects of laminar separation may appear and 
some pressure drag might occur. The small differences 
in drag between the, envelope polar for the new airfoils 
and for the NACA 2l+69-5ij. airfoil are probably a result 
of this phenomenon. Actually tiie point of_ maximum nega- 
tive pressure for the new airfoils is considerably 
farther back than the corresponding point for the NACA 
2ii.09-5i4. airfoil but, if laminar separation occurs early, 
nearly equal drag coefficients might be expected. 

At high speeds (M = O.75, fig. 20), the region Tor 
which the NACA l6-series airfoils were" desig^ied, the 
superiority of the new airfoils is clea'r.' ~The"' "earlier 
onset of the compressibility ef fects ,f or the older air- 
foils leads to early drag increases and lowered maximum 
lift coefficients. At speeds above M = 0.75 the use of 
the 6lder sections appears unwarranted for any purpose. 

Critical speed .- The variation of the critical 
speed with lift coefficient and with thickness is given 
in figures 21 and 22, respectively. " These curves indicate 
that critical speeds exceeding the theoretical values 
were attained in the tests. In the choice of the test 
critical speeds, the values were selected on the basis 
of earlier experience that indicated some rise in drag 
before large flow disturbances occurred. If these speeds- 
were chosen as the highest values reached before any 
appreciable drag; increment occurred, the agreement with 
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the theoretical .curves would be very gao'd. For ooniparlson 
the critical speed of the 3c8 aii-roll is plotted. in £a.g-- 
ure 21". The difference between .the new and the older 
airfoils is greater than shown by the curves because the 
5C3 is 8 percent thick, or 1 percent of 'the chord thinner 
than the airfoils pf the NACA I6-609 series. 

Minimum drag . - Coefficients of minimum drag plotted 
ageilnst Reynolds number are given in figures ,25 and 2J4.. 
^he lowest drag coefficient was obtained for the NACA 
16-106 airfoil; this coefficient is approximately 0,0026 
at low speeds and increases ttr-approximately O.OO32 
immediately below the critical gpeed. Cf the 9-percent- 
thick airfoil series designed to operate at various lift 
coefficients^ i?he NACA I-6-IO9 airfoil appears to have 
the lowest drag. This result is contrary to expectation 
because the symmetrical or basic form of the NACA 16-OO9 
would normally have the lov/est minimum drag ooeff iciient. 
The "difference may be due to some irregularity of the 
airfoil surface. 

The comparison of the. minimum drag coefficients 
for the 3C8 and the NACA 2i|.09-3i^. and 16-209 airfoils is 
shown in figure 25 . The high critical speed -for the NACA 
16-209 airfoil is apparent.. The comparison as given- 
directly by figure 25 is a little misleading because of 
the smaller thickness ratio for the JCS airfoil. For 
equal thickness ratios, the differences between the 
C-series and the NACA l6-series airfoils will be greater 
than shown. 

Use of the data . - The envelope polars that may be 
drawn for the NACA I6 -series airfoils represent a "new 
and much lower drag as well as higher critical speed "■ 
attainable for the design of propeller-blade sections. ' 
Even though the angle-pf -attack range Is less than for 
the older sections, there will be numerous designs for' 
which sufficient angle-of -attack range is given by the 
new sections, ^ov high-speed, high-altitude aircraft, 
the advantages of the low drag and hig^i critical speed 
are of paramount Importance and, in these designs, rational 
choice of section is of increasing" Importance . In many 
designs the diameter is fixed by considerations other ■ 
than p3r»opeller efficiency. Thus the induced losses are 
fixed and propellers of hlgjiest efficiency can be developed 
only by operating and designing- the blade sections to 
operate on the envelope polfsr^s. Another important con- 
sideration In using new blade sections to achieve highest 
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efficiency concerns the adaptation of the sections to 
older propeller designs. Optimum efficiency cannot be 
achieved by simply substituting the new sections for the 
old on a givBn design. The use of better blade . sfeetions 
permits the use of larger diameter and necessitates some 
plan-form changes., .All these factors should be .considered 
in a design for bes't" eTf iciency with the new blade" sections. 



CONCLUSION 



By a new approach to airfoil design based upon 
findings of fundamental flow studies, a new series of 
airfoils, the NACA I6 series, have been develoged which 
have increased critical Mach number and at low speeds 
reduced drag. ... .. . - ,"- _ Jr-:~-^■^-.•-~^- 
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TABI^I I 



CAMHSR-LINE ORDINATES FOR NACA I6- AND 07-S3RIES 
AIRFOILS V/ITH = -L.p ^ 

[All values measured in percent chord 
from chord Ixn^ 



Station 



0 

1-25 

^ 5. Oft 

7.5^ 
,10.00 

.15.00 

20^00 

.25.00 

; ^Q.oo 
14.0-00 
50-00 

60 -00 
70 .00 

,90.00 
95-06 
100.60 



Ordinate 



0 

.930 
1.580 

2.120 
2. 587 

5.556 

5,516 

4,861 
5.982 

2.5G7 
1.5S0 



Slope 



^3gT. 
.2^155 

.lSfe95 

,m66 

.15601+ 
. 110^2 

.087U3 

.0671+5 
.03227 

0 

- .03227 

- .067i+5 
^ .1105 2 

- .17W 

- .25i|-52 

- .6225I4. 
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THICKNESS ORDI NATES FOR AIRFOILS WITH 



THICKNESS 9 PERCENT OF CHORI)-'- 

[All ^''values measured in percent chord from 
I and perpendicular to camber llnej 



Station 


Ordinate s 


NACA 16-serles 

airfoils 


NACA 07-009 

airfoil 


0 

2*5 .V , 

5 s 
7.5 r.H 
.10 • 

20 , 

50 

ko 

.50 

.6p . ■ 

.io ■:■ 

50 

9:5 

100 


b*1.55li.' ■'■^'^ov 
■ fl.BSi i"'^''' 

•'2.593-" 

.3.1Ql-_. v.: -'ill 
;5.U98"^' 


0 

1.25 
1.67 
2.19 
2.58 
2.90 

5.41 
5.79 
4.07 

4.27 

4.46' 

4.50 
4.37 

k.OQ 

5.54 
1.91 
1,00 

.09 


.4.065^ v"{2i 
4.591^- t.--? 

. -'Il.5.oo.i< 

■4.576V- ^■■"'-^^ 
5.952 

5.149-' --i-'^' 

1.886V-" 
.1.061 

.:.090 


Slope of radius «4^^3-4 f 
througji .^rii^ of ahord *= 0ibaa>40L ' 1 


, L.E. radius of 

NACA 16 -aeries airfoils = 0.5.96 ( t/O.09 )^ 



B'or other thicknesses (t, in perpent) 
multiply ordinates for NACA l6-series air- 
foils by t/0.09.: 
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Figure 1'.- Basic airfoils and theoretical pressure 
distributions . 
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Figure 3.- Comparison of minimum drag obtained by 
various methods. 
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Figure 2.- Profiles for airfoils having high critical speeds 
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Figure 5.- Aerodynamic cliar&cteristios of BTAGA 16-009-serieB 
airfoils. M « 0.45. 
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(a) Polar plots. 



(b) Lift aad moment data. 



Figure 5.- Aerodynamic cliaracteristios of NAGA 16-009-serie6 
airfoils. M « 0.45. 
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Figure 7.- Aerodynamic characteristics of NAOA 16-009-series 
airfoils. M « 0,70. 
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Figure 9.- Aerodynamic charaot eristic a of NAOA IS-BOO-serles 
airfoils. VL a 0.50. 
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Figure 10.- Aerodynamic characteristics of NAOA 16-500-series 
airfoils. M = 0.45. 
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Figure 11.- Aerodynamio cliaracteristics of NAOA 16-500-Berie8 
airfoils. M » 0.60. 
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Figxire 13.- Aerodynamic characterietioe of NACA 16-500-8eries 
airfoils. U = 0.70. 
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Figure 13.- Aerodyneanio charaoteristlce of HACA IS-SOO-series 
airfoils. . M » 0,75. 
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Figure 14.- Aerodynamic characteristic 8 of NAOA IS-OOO-series 
airfoils. U a 0.20. 
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Fipnire 15.- Aerodynamic oharacteristiCB of NAOA IS-OOO-series 
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Figure 16.- Aerodynamic characteristics of NACA IS-OOO-series 
airfoils. H » 0«60, 
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Figure 17.- Aerodynamic oharaoteristics of NACA 16-000-serie 
airfoils. M « 0.70, 
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Figure 18.- Aerodynamic characteristics of NAOA IS-OOO-eeries 
airfoils. M = 0,75. 
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Figure 19.- Comparieon of airfoils. M » 0«45« 
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Figure 20.- Comparison of airfoils. 21 s 0.73. 
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Flgiore 21.- Variation of airfoil critical Uach number with, 
lift coefficient. 
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Figure 22.- Variation of airfoil critical Mach number with 
thickness-chord ratio. 
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Figure 25.- Comparison of airfoil minimxua drag coefficients 



